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SYSTOLIC  ARRAY  AOAPTIVE  BEAMFORKING 


INTRODUCTION 


Optimum  algorithms  for  the  space-time  processing  requirements  of  a 
discrete  sensor  array  have  existed  1  i  one  form  or  another  for  almost 
twenty-five  years  [Ref.  1].  The  computational  requirements  for 
Implementation  have  Inhibited  the  widespread  application  of  these  techniques 
to  broadband  arrays.  To  date,  the  hardware  realization  of  the  Intensive 
linear  algebra  operations  required  for  any  type  of  modem  slgna1  processing 
function  has  not  been  amenable  to  cost  effective  solutions.  In  this  regard, 
the  most  promising  new  development  In  modern  signal  processor  design  has 
been  the  result  of  rapid  progress  In  very  large  scale  Integration  (VLSI) 
technology.  Large  scale  Integration  has  allowed  the  fabrication  of  special 
purpose  components  which  has  provided  the  Impetus  for  the  evolution  of  very 
powerful  special  purpose  computing  architectures  for  linear  algebra 
Intensive  signal  processing  requirements  [Ref.  2].  While  this  field 
represents  an  area  of  current  Intensive  research,  some  concepts,  such  as  t.he 
systolic  computing  cellular  array,  have  already  produced  significant 
developments  [Ref.  3]. 

The  systolic  array  and  related  wavefront  processor  are  specifically 
designed  to  exploit  the  unique  regularities  of  a  particular  linear  algebra 
operation  [Ref.  4].  In  particular  a  characteristic  of  many  matrix  algebra 
operations  Is  the  requirement  for  data  communication  between  only  nearest 
neighbor  arithmetic  cells  In  a  properly  designed  array  of  such  cells.  This 
basic  principle  of  simplification  in  conjunction  with  the  relatively  low 
cost  of  VLSI  arithmetic  cell  components  makes  It  feasible  to  design 
hardwired  systolic  array  algorithms  with  essentially  no  internal  control, 
minimal  memory,  and  maximal  parallelism.  As  a  case  In  point,  this  report 
discusses  the  broadband  minimum  variance  distortionless  response  (MVDR) 
beamforming  algorithm  that  consists  of  the  following  three  matrix 
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operations:  Cholesky  factorization  of  an  estimated  cross-spectral  density 
matrix  (CSOM),  solution  of  a  least-squares  (filter)  problem  after  each  rank 
one  update  of  the  CSOM,  and  an  N-channel  matrix  filter  operation. 

First,  the  theory  and  direct  element  level  Implementation  of  the 
adaptive  broadband  MVDR  beamformer  is  presented,  then  the  systolic  array 
Implementation  Is  described.  This  Is  followed  by  a  discussion  of  the 
theoretical  performance  predictions  for  an  MVDR  process.  Finally,  some 
Implications  of  systolic  array  architectures  with  respect  to  variations  of 
the  MVDR  algorithm  for  high  resolution  space-time  processing  In  very  large 
arrays  are  considered. 


THEORY  AND  DIRECT  IMPLEMENTATION 


A  convenient  discrete  frequency  domain  representation  of  the  broadband 

data  from  an  N-sensor  array  Is  given  by  the  vector  x.  with  transpose 

T  * 

x^  specified  by 


5k  *  txlk<wl*  JWwl)xlk(w2)  xNk(w2>  “•  xlk(“N)  xNk(wM^ 


r  T  T  T  . 

^lk  x2k  •••  ’ 


(1) 


wnere 


x1k<“»>  *  jr  j 


-jut 


m 


x^tje  dt  , 


(2) 


tk  '  2 


with  x^(t)  the  output  time-domain  waveform  of  the  1-th  sensor  and  uffl  the 
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m-th  discrete  radial  frequency 


»  -  2vm/T  0  <  m  <  W-l . 

f!i  —  - 


(3) 


Therefore,  the  elements  of  the  vector  are  seen  to  be  the  discrete 
Courier  coefficients  of  x^(t)  over*  the  Interval  (t^  -  T/2,  tk  T/2). 

When  the  observation  time,  T,  Is  large  with  respect  to  the  Inverse  bandwidth 
of  x^(t),  the  Fourier  coefficients  between  frequencies  become 
uncorrelated.  We  shall  assume  that  this  condition  Is  satisfied  and  that  the 
waveforms  x^(t),  1  <  1  <  *►,  while  not  necessarily  fiausslan  random 
processes,  are  zero-mean  and  wide  sense  stationary. 

Consistent  with  the  above  waveform  Assumptions,  the  covariance  matrix 
R  *  R(w)  at  frequency  «  ,  hereafter  referred  to  as  CSDH,  '.s  defined 
with  Ij-th  element  Et*tk(“w)xjk(«m)]*  The  covariance  matrix  for 
the  Fourier  coefficient  vector  Is  an  MN-by-MN  block  diagonal  matrix 
with  the  CSDM  R^  as  the  m-th  N-by-N  main  block  diagonal  element.  The  MVDR 
beamformer  filter  at  time  tk  Is  now  defined  In  terms  of  the  complex 
conjugate  transpose  (Indicated  by  suoe'scrlpt  H)  as 


-k  “  twlk(wl) 


wNk<wl )wik(“2)  wNk*w2'  *'•  tolk(wM>  •••  wNk<wH^ 


-  c«; 


H  H 
kl  -k2 


n  *1 

*kMJ  * 


(4) 


The  total  broadband  output  power  of  the  MVDR  beamformer  is  defined  as  the 
variance 


2  CM  H  I21 

°  *  £t|*k  J 


H  H. 

«k  E[xk  xkj  wk 
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M 

"  ^  ^mk  Rm  ^mk  ‘ 

m-1 

The  requirement  for  a  distortionless  spectral  response  to  e  signal  with  the 
steering  vector  d  «t  frequency  Is  equivalent  to  the  constraint 


Refwj^  d^]  -  l  and  Iwfw^  d^]  ■  0,  1  <  ■  <  N  (6) 

on  the  MVDR  filter  weight  vector  w.  .  The  steering  vector 

n-th  element  exp  (j*»r  )  where  r  Is  the  relative  time  delay  to  tne 

m  op  np 

n-th  sensor  for  a  signal  from  the  p-th  direction.  Using  the  method  of 

2 

undetermined  Lagrange  multipliers,  a  of  Eq.  ($)  Is  minimized  subject  to 
the  constra4nts  of  Eq.  (6)  If  the  criterion  function 


•m  has  the 


M 

*  * IJ  [£  %  **  ♦  2k»*  V  -1  >  *  *h.  < «*t«2k  v 1  (,) 


Is  minimized  with  respect  to  the  real  and  Imaginary  parts  of  v» 
simultaneously.  The  solution  to  this  problem  Is  known  to  be 


«mpk 


rfl/d" 

-mp  nap 


R« 

m  -mp 


1  <  m  <  M 


(8) 


where  the  beam  Index  p  Indicates  that  a  unique  beam/f liter  vector  w  Is 

“TOP  It 

required  for  each  of  the  B  (1  <  p  <  B)  beam  steering  signal  directions.  The 
corresponding  minimum  total  beamformer  broadband  output  power  (variance) 
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obtained  using  Eg.  (8)  In  Eq.  (5)  is 


bp  (-t«P  R*  %)"1,  (0) 

i»-l 

Therefore,  the  estimated  signa’  autopower  spectral  density  at  frequency  Is 

'sUpC  v*1-  no, 

In  pra-Mce,  the  CSOM  matrix  cannot  be  estimated  exactly  by  time 
averaging  * jeause  the  random  process  *n(t)  Is  never  truly  stationary 
ana/or  ergodlc.  As  a  result,  the  available  averaging  time  Is  limited. 
Accordingly,  one  approach  to  the  time-varying  adaptive  estimation  of  Rr  at 
time  t^  Is  to  compute  the  exponentially  time  averaged  estimator  of  the 
CSDM  R  at  time  t.  as 

fl  K 


%*  -  «Vw  *-2  •  <"> 

where  *  Is  a  smoothing  factor  (0  <v<l)  that  Implements  the  exponentially 
weighted  time  averaging  operation.  Eq.  (11)  Is  a  rank  one  update  time 
averaging  operation  to  the  CSDM  at  frequency  w  and  requires  approximately 

7  I" 

MR n  complex  multlply-addltion  operations  and  an  equivalent  memory  size. 

The  estimator  R  Is  used  In  Eqs.  (8),  (9),  and  (10)  to  compute  the  MVDR 

INK 

filter  weights,  broadband  beamforme*'  output  power,  and  estimated  signal 
spectrum,  respectively.  A  straightforward  Inversion  of  R^  requires  on 
the  order  of  Mh3  operations.  Approximately  BMN2  additional  operations 
are  required  to  form  the  beam  filter  vectors  In  Eq.  (8)  and  beampower 
estimates  In  Eq.  (9)  for  a  system  ulth  B  beams.  The  average  number  of 
operations  for  the  direct  Implementation  of  the  broadband  MVDR  beamformer 
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is,  therefore,  on  tho  order  of 


0d  -  M((d2/2)  ♦  8N2  ♦  (8N2  ♦  N3)/K) 
«  MN2  (B  ♦  (B  +  N)/K) 


(12) 


per  beamformer  update  cycle.  The  parameter  K  Is  the  number  of  update  cycles 
of  the  estimated  CSDM  R_  between  computations  of  the  updated  filter 

III 

weights  and  beam  output  power  estimates. 


SYSTOLIC  IMPLEMENTATION 


In  terms  of  systolic  computing  array  functions.  Implementation  o,!  the 
MVDR  beamformer  does  not  substantially  reduce  the  computational  burden 
Indicated  by  Eq.  (12).  Instead,  systolic  computing  exposes  the  regularity 
of  a  particular  function  In  such  a  way  as  to  make  that  function  amenable  to 
the  maximum  of  parallel  computation  and  a  minimum  of  both  control  and  data 
transfer.  In  the  following,  the  primary  systolic  computing  elements  are 
developed. 

First,  consider  the  requirement  to  compute  the  estimated  output  power 
for  the  m-th  frequency,  p-th  beam  steering  direction,  and  k-th  update  cycle 
from  the  beam  output  power  estimate 


bmpk  *  1/gmpk 

where 

H 

^mpk  *  ^mp  Rmk  ^mp* 


(13) 
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In  the  systolic  Implementation  to  be  described  the  Inverse  CSDM  rJ  Is 

INK 

never  formed  explicitly.  Rather  the  upper  triangular  Cholesky  factor 
of  R^,  defined  by 

UK 


no 


Is  obtained  directly  (as  explained  later).  The  diagonal  elements  of  the 
Cholesxy  factor  are  real.  Substituting  Eq.  (14)  Into  (13)  yields 


*mpk 


«C>'V  vH  “O'1  v 


•■mpk  nepk 


I  ^mpk  I 


(15) 


H  -1 

where  only  the  complex  vector  v  .  -  (U  . )  d _  must  be  computed  at  a 

0  iripK  mK  "trip 

cost  of  HN  /2  operations.  To  form  the  beam  output  from  the  beam  filter 
vector  of  F.q.  (8),  i.e., 


'mpk 


%npk  ^mk 


^mp  Rmk  *mk/9mpk' 


again  use  the  Cholesky  factorization  of  Eq.  (14)  to  obtain 


w  ■  «&>''  vH  «ul<r'  %''V 


H  H  2 

\pk  ^mpk'l^mpkl  ' 


(16) 
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The  Important  point  with  respect  to  Eq.  (16)  Is  that  the  computation  of  the 
vector 


^mpx  *  ^Umk^  ^*nk 


Is  of  the  exact  form  as  that  for 


,„H  v-1  , 

^mpk  =  ^Umk^  - 


In  Eq.  (15).  Thus,  the  same  circuitry  can  be  used  for  each  computation. 

A  realization  of  the  MVDR  beamformer  as  described  above  requires  the 
Implementation  of  two  functions  with  systolic  arrays:  a  rank  one  update  of 
the  CSDM  Rffl)(  In  terms  of  Its  Cholesky  factor  Umfc  In  Eq.  (14)  and  the 
subsequent  solution  of  a  set  cf  N  linear  equations  of  the  form 


u!_ .  h  =  b 

mi.  -  — 


as  required  by  Eqs.  (17)  and  (18). 


Rank  One  Cholesky  Factor  Update  T 5 . 6 1 :  Let  the  N-dlmenslonal  row  vector 


.1/2  H 
I  -  (1  -  v)  2L. 


'  [Z1  z2  zN! 
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and  the  upper  triangular  H-by-N  matrix 


U  -  v 


1/2 


m,k-1 


(22) 


be  defined.  Now  form  the  (N  +  1)-by-N  data  augmented  Cholesky  matrix 


(23) 


which  can  be  reduced  to  upper  triangular  form  by  a  sequence  of  N  orthogonal 
plane  rotations  Indicated  by  the  orthogonal  premultlplleir  matrix  P  as 


Z1  Z2 


U11  U12  ‘ 


0  U22  * 


•  Z- 


.  u 


IN 


.  u 


2N 


0  0  ...  u 


NN 


U11  u12  U1N 

0  u22  U2N 

o  o  •••  «NN 

0  0  0 


(24a) 


(24b) 
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where  P  *  Pu  P,.  ,  ...  P, ,  P„  P*  I  and  o  is  an  N-vector  of  all 
n  n-i  inn  - 

zeroes.  Thus, 


BHPHPB  -  BHB 


-  UHU  +  zHz 


<k-l  Vk-1  ♦  < 1  -  •>im 


«  u 


H 

mk 


(25a) 

(25b) 

(25c) 

(25d) 

(25e) 

(25f) 


and  U  Is  the  rank  one  updated  Cholesky  factor  of  the  estimated  CSDM  Rmk. 

If  1  <  1  <  N  ts  real,  then  u^,  1  <  1  <  N  Is  also  real  when  the 
boundary  cell  process  definition  of  figure  i  Is  used. 

The  systolic  array  realization  of  the  rank  one  Cholesky  factor  update 
requires  both  a  boundary  cell  and  an  Internal  cell  as  Illustrated  In  figure 
1.  The  specific  linear  systolic  array  for  reduction  of  Eq.  (23)  to  upper 
triangular  form  Is  given  In  figure  2.  It  Is  the  boundary  cell  that  must 
compute  the  correct  plane  rotation  angle  complex  cosine  and  sine  values 
(c,s)  and  the  Internal  cells  which  propagate  tnat  rotation  down  the  n-th  row 
of  the  matrix  Pn_j  Pn2  ..  P1  U. 


Linear  Equation  Back  Solve  (71:  The  solution  of  Eq.  (19)  is  performed 
a.  a  sequential  backsubstitution  of  the  elements  of  h  as  they  are  obtained 
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t 


* 


i 


input  . -  '  >  process  — —  -->»  output 

(plane  rotation) 


Figure  1A.  Boundary  Ceil  Input-Process -Output  Diagram 
for  Cholesky  Factor  Update  (real  y) 


Delay: 


9 

input 


(delay) 

process 


output 


Compute: 


input 


process 


- output 


Figure  IB.  Internal  Cell  Input -Process -Output  Diagram 
for  both  the  Delay  and  Compute  Function  Required 
for  Cholesky  Factor  Update 


Figure  1.  Boundary  Cell  and  Internal  Cell  Input  Process -Output 
Diagram  for  the  Cholesky  Factor  Update 
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1  ITS 


H2 


u22  *  U13 


u23  W14 


133  24 


(>t^Xp  ■  =0 

•  Cy  •  •  •  •  • 


I 

*11  ‘  Z2 

u12  '  Z3 

J22  *  u13  * 

U23  *  U14 


•  •  •  • 


'NN 


Figure  2.  Linear  Systolic  Array  for  a  Rank  One  Cholesky  Factor  Update 
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4 


in  sequence.  Eq.  (19)  Is  of  the  form 


un  0  0  ...  0 

V 

hi 

u12u22°  ...  0 

h2 

- 

b2 

“13  u23  “33  ••  ' 

«  •  •  « 

• 

h3 

b3 

• 

•  •  •  •  ♦ 

“in  U2N  U3N  *'•  UNN. 

k 

- 

> 

Thus 

"v '  Vun 

h2  *  (hp  ■  ui2^1  ^u22 
•h3  =  *b3  “  U13hl  •“  U23*V/U33 

hN  =  (bn  "  UlNhl  **•  UN-l,NhN-l)/uNN  ’ 


(26) 


(27) 


which  Is  realized  with  an  Implementation  of  the  linear  systolic  array  using 
the  cells  Illustrated  In  figure  3.  The  corresponding  systolic  array  is 
shown  In  figure  4. 

It  Is  noted  that  systolic  arrays  for  the  rank  one  Cholesky  update 
(figure  2)  and  linear  equation  backsolve  solution  (figure  4)  are  linear 
arrays  with  minor  differences  In  input-output  structure.  As  discussed 
above,  utilization  of  the  boundary  and  internal  cells  is  fifty  percent.  One 
hundred  percent  utilization  is  achieved  with  the  interleaving  of  two 
problems  delayed  by  one  cycle.  Moreover,  it  is  possible  using  a 
nonalgebraic,  i.e.,  geometric  approach,  to  partition  large  (N)  problems  into 
segments  wherein  a  smaller  systolic  array  of  size  less  than  N  internal  cells 
can  be  time  multiplexed. 
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Input  process  . output 


Figure  3A.  Boundary  Cell  Input-Process -Output  Diagram 
for  Linear  Equation  Backsubstltution 


Figure  3B.  Internal  Cell  Input-Process-Output  Diagram 
for  Linear  Equation  Backsubstltutlon 


Figure  3.  Linear  Equation  Backsubstitution 
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Figure  4.  Systolic  Array  for  Solution  of  Linear  Equations 
by  Backsubstltutlon  for  a  4-Element  Partition 
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PERFORMANCE  OF  AN  MVOR  BEAMFORMER 


The  CSOM  for  an  array  that  receives  both  a  signal  with  autopower 

2 

spectral  density  (ASO)o  and  an  Interference,  which  Is  uncorrelated 

5  2 

with  the  signal,  and  an  ASO  at  frequency  f  Is  given  by 


*  -  A  4,  A  *  A  4, 4?  ♦  A  i.  <»> 

2 

In  Eq.  (28),  aQ  Is  the  spatially  uncorrelated  noise  ASD  at  frequency  w, 

I  Is  an  N-by-N  Identity  matrix,  and  and  are  the  steering  vectors 
for  the  signal  and  Interference,  respectively.  The  interference  Is  assumed 
to  be  In  the  farfleld  of  the  array,  l.e..  It  Is  a  Point  INterference  (PIN) 
with  respect  to  angle  of  arrival. 

The  MVOR  filter  vector  for  beamformlng  In  the  signal  direction  at 
frequency  w  has  been  shown  to  be 


w  -  R'1  d$/d”  R"1  ds  .  (29) 

For  a  conventional  time  delay-and-sum  beamformer  (CBF)  the  vector  d^ 

Is  used.  The  gain  In  signal  to  the  total  of  Interference  and  noise  power  at 
the  beam  output  relative  to  a  single  sensor  is  given  by 


gmv  *  +  +  “lu  (30) 

for  the  MVOR  beamformer  and 

Gc  =  >  o*I]ds)/(o^/[^  4-  oj])  (31) 
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for  the  CBF.  The  array  gain  Improvement  (A6I)  of  the  NVDR  relative  to  the 
CBF  Is 


Aft!  m  ft  /ft 

noi  >V'°C 


-  1  ♦ 


V-r 


(32) 


where 

r  -  No2/*2  (33) 

and  represents  the  Interference  to  uncorrelated  noise  variance  ratio  at  the 
output  of  a  CBF  steered  with  Its  main  response  axis  (MRA)  directly  at  the 
Interference.  The  parameter  I  (o  <  |  <  1)  Is  defined  by  the  expression 


1  -  uJitlW.  (34) 

Figure  L;  gives  the  AG!  as  a  function  of  r  for  four  different  values  of  1; 
namely,  I  -  1/2  (-3  dB),  I  -  2/100  (-17  dB),  I  -  5/1000  (-23  dB)  and  1  * 
1/1000  (-30  dB).  For  either  high  Interference  levels  and/or  low  noise 
levels,  significant  gains  tf  greater  than  10  dB  are  achievable  with  MVDR 
beamtormlng  over  a  C'lF  At  high  noise  levels  the  Increased  urcorrelated 
noise  begins  to  mask  the  PIN  as  the  major  degrading  factor  ard  useful  AGII 
Is  only  achieved  wnen  the  PIN  Is  Inside  the  mainlobe  of  the  beam  receiving 
the  signal  (-3  dB  re  MRA) . 


AOAPTIVE  BEAMFORMING  FUR  VERY  LARGE  ARRAYS 

The  fundamental  parameter  whl-ch  determines  the  signal  to  total 
background  noise  variance  ratio  gain  for  an  array  is  the  number  of  sensors 
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« 


Figure  5.  Minimum  Variance  Optimum  Beamformer  AGI  Relative 
tc  a  CBF  as  a  Function  of  the  Interference-to-Nolse 
Ratio  for  a  PIN  Measured  at  the  Output  of  a 
CBF  Steered  Directly  at  the  Interference 
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fi.  SlMply  stated,  If  an  array  cannot  provide  a  certain  level  of  required 
conventional  beamformer  a  ray  gain  with  only  spatially  uncorrelated  white 
noise  present,  then  adaptive  beamforming  will  not  alter  this  fact  and  the 
only  option  to  make  the  number  of  sensors  (N)  In  the  array  larger.  To 
Illustrate  this  point  It  Is  observed  that 


®NVl  2  *®c|?  .  (J5> 

^  ■  0  ■  0 

-  N  . 

Thus,  from  the  array  gain  perspective,  Cb.-  Is  optima*  when  there  Is  no 
Interference  present  and  the  only  recourse  Is  to  build  larger  arrays,  l.e., 
arrays  with  itr>re  sensors, 

Given  the  specific  spatially  uncorrelr.ted  white  noise  array  gain,  N,  of 
Eq.  (35)  for  an  array  of  N  sensor  elements,  an  Implementation  of  the  MVDR 
beamformer  described  previously  requires  the  update  of  an  N  dimensional 
Cholesky  factor  and  the  solution  of  a  correspondingly  large  set  of  linear 
equations.  Thus,  the  number  of  sensor  elements  N  equal  to  the  array  white 
noise  gain  determines  the  size  of  the  element  space  MVDR  system.  For  arrays 
with  a  large  number  of  sensor  elements  N,  the  computational  requirements  can 

3 

become  prohibitive  given  that  the  computing  burden  Is  proportional  N  as 
specified  by  Eq.  (12). 

In  a  dynamic  situation,  where  the  angle  of  arrival  for  a  particular 
Interference  Is  changing  with  time,  the  effective  averaging  time  for 
estimating  the  Interelement  CSOM  Cholesky  factor  Is  limited  by  a  temporal 
ststlonarlty  assumption.  Thus,  the  variance  of  the  elements  In  the  CSDM 
estimator  of  Eq.  (11),  which  Is  Inversely  proportional  to  averaging  time, 
has  a  lower  bound  determined  by  the  finite  averaging  time.  Specifically,  If 

M  Is  the  effective  number  of  statistically  Independent  sample  vectors,  x  , 

m 

which  are  exponentially  averaged  to  produce  the  estimated  CSDM  of  Eq.  (11), 
then  the  variance  on  the  MVDR  beam  output  power  estimator  detection 
statistic  of  Eq.  (10)  Is  Inversely  proportional  to  M  -  N  +  1  where  It  Is 
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assumed  that  M  >  N  [Refs.  8,  9].  Thus,  as  the  number  of  elements  (N)  In  the 
array  Increases,  It  Is  necessary  to  Increase  the  CSOM  estimator  averaging 
time  as  determined  by  K  proportionately  to  maintain  the  same  beam  output 
power  estimator  variance. 

eventually,  for  the  element  space  MVDR  process,  the  size  of  the  array  N 
Is  limited  by  the  time  statlonarlty  constraint  which  Is,  In  turn,  determined 
by  the  Interference  position  rate  of  change  with  respect  to  time. 

A  natural  way  to  avoid  the  temporal  statlonarlty  limitation  on  the  white 
noise  array  gain  N  discussed  above  Is  to  perform  the  systolic  MVDR  process 
In  a  domain  other  than  the  N-dlmenslonal  element  space.  If  this  new  domain 
has  a  lower  dimensionality,  both  the  systolic  engine  computational  and 
memory  size  complexity  and  the  effective  time  averaging  requirement 
constraints  are  reduced  accordingly.  References  [10]  and  [11]  suggest  the 
liNplefiKi'.tatlo.i  of  the  MVDR  «ilgor1*hm  In  e  so-called  beam  space.  In  beam 
space,  only  spatially  orthogonal  CBF  beams  that  are  steered  contiguous  to  a 
selected  reference  beam  are  used  as  inputs  to  an  MVDR  beam  Interpolation 
algorithm.  Clearly,  the  question  of  selecting  the  appropriate  number  of  and 
location  for  orthogonal  beams  Is  not  straightforward.  At  a  minimum,  this 
selection  Is  a  frequency  dependent  process  due  to  the  variation  of  beam 
overlap  caused  by  the  Increase  of  beamwldth  with  a  decrease  In  frequency. 

In  addition,  the  number  of  Independent  beams  must  be  made  large  enough  to 
provide  a  sufficient  number  of  degrees  of  freedom  for  near  optimum 
performance  In  a  multiple  Interference  condition. 

As  a  practical  matter,  even  the  formation  of  a  conventional  time 
delay-and-sum  beam  for  a  very  large  array  Is  a  difficult  implementation 
Issue.  Usually  partial  aperture,  l.e.,  subarray,  beams  are  formed  as  a 
first  step  In  the  formation  of  a  full  beam  from  a  large  array.  An 
alternative  to  the  beam  space  approach  for  dimensionality  reduction  in  very 
large  arrays  Is  referred  to  herein  as  the  subarray  (SA)  space  formulation 
[Refs.  12,  13].  In  this  approach,  subapertures  of  contiguous  elements  in 
the  large  array  are  prebeamformed  using  simple  time  delay-and-sum  and  fixed 
spatial  windowing  techniques.  This  partitioned  subarray  beamforming  (SBF) 
can  be  envisioned  as  creating  a  secondary  array  of  spatially  directional 
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elements  which,  In  turn,  ere  processed  with  an  N/P-dlmenslon  NVDR  beamformer 
In  cascade  with  the  subaperture  beamformers  discussed  above. 

If  each  of  the  partitioned  subaT&ys  (N/P),  as  Illustrated  In  figure  6, 

consists  of  P  contiguous  sensor  elements,  then  the  MVDR  process  Is  of 

dimension  N/P.  However,  as  with  the  beam  space  approach,  more  than  one 

small  (dimension  N/P)  CSDM  Cholesky  factor  needs  to  be  estimated  at  each 

frequency.  This  Is  In  contrast  with  the  element  space  NVDR  where  a  single 

very  large  (dimension  N)  CSON  Cholesky  factor  Is  estimated.  This  Is  because 

■  each  SBF  can  form  approximately  P  spatially  Independent  beams  which  resolve 

substantially  nonoverlapping  segments  of  solid  angle.  Thus,  for  the 

formation  of  a  particular  CSON  matrix  estimator,  those  SBF  outputs  steered 

at  the  same  angle  should  be  selected.  This  SA  space  requirement  would 

constitute  a  need  for  approximately  P  NVDR  parallel  processes  each  of 

dimension  N/P  as  opposed  to  one  NVDR  process  of  dimension  N  required  for  the 

element  space  formulation.  It  Is  noted  that  from  the  computational 

requirement  standpoint  the  SA  space  burden  Is  proportional  to  B,..  ■ 

2  3 

(N/P)  N  as  contrasted  to  BE  *  N  for  element  space.  The  actual 

burden  Is  proportional  to  ((N/P)2  P  ♦  (N/P)2  N]  (N/P)2  N  for  N  » 

P.  The  Cholesky  factor  update  burden  Is  (N/P)2  P  and  the  backsubstltutlon 
2 

burden  Is  (N/P)  N.  Thus,  the  computational  load  and  memory  size 
reductions  can  be  enormous  when  the  SA  space  Is  adopted.  Moreover,  the 
restriction  on  the  effective  averaging  time  N  Imposed  by  the  array  size  N 
becomes  N  -  (N/P)  ♦  1  >  T,  where  T  Is  a  threshold  set  by  the  desired 
variance  of  the  beam  output  power  detection  statistic.  It  follows  that 
averaging  time  can  be  reduced  In  a  SA  space  formulation  to  accommodate  the 
*  spatial  dynamics  of  the  Interference  with  essentially  no  loss  of  performance. 

The  beam  and  SA  space  MVDR  array  gain  performance  Is  obtained  by 
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NELEMENT 
ARRAY 


CHANNEL 
* 

BEAMFORMER 


FIRST  STAGE 
BEAMFORMERS 


SECOND  STAGE 
BEAMFORMERS 


Figure  6.  Subarray  Space  Partitioning  and  Cascaded  Beamforming  of  a 
Very  Large  N-Element  Array  to  Reduce  N/P-Input  MVDR  Beamformer 
Complexity  and  Convergence  Time  Restrictions.  (Each  of  the  (N/P) 
subarrays  consists  of  P  elements  with  conventional  beamforming.) 
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Introducing  the  array  element  data  preprocessing  matrix 


» 


Dx 


P- 


L/2  . . 


*sl  ft 


ft  *s2 


^0 


.  for  beamspace 


(36a) 


0 


fts,N/p) 


for  SA  space. 


(36b) 


In  Eq.  (36a),  d^  is  a  C8F  N-dlmenslonal  steering  vector  corresponding  to  a 
beam  space  patch  of  size  L+1  which  Is  defined  as  being  centered  at  a  point, 
e,  specified  by  a  reference  beam  which  has  steering  vector  d  «=  d  (e). 

U  U  y 

Ideally,  these  steering  vectors  are  assumed  to  be  orthogonal,  l.e.,  d" 
dj  =  N4ij-  For  the  SA  space  MVOR  process,  the  P-dlmenslonal  vector 
ftsk  corresponds  to  the  P-element  subvector  of  the  N-dlmenslonal  steering 
vector 


j 


-S.N/P 


(37) 


which  would  be  required  to  electronically  steer  the  entire  array  subarray  by 
subarray  at  a  single  point  that  Is  Implicit  In  the  steering  vector  dQ. 

For  the  beam  space  formulation  0  Is  an  N-by-(L-t-l)  dimensional  matrix  and  for 
the  subarray  counterpart  0  Is  of  dimension  N-by-(N/P).  If  (L+1)  =  N/P  for 
these  two  suboptimum  MVDR  processes  and  the  same  AGI  performance  results, 
then  the  two  approaches  would  require  the  same  hardware  Implementation  with 
identical  performance  except  that  the  beamspace  process  requires  full 


23 


TR  7981 


conventional  beams  to  be  formed  Instead  of  only  partial  beams. 

To  establish  the  MVDR  AGI  for  the  two  suboptimum  procedures,  the  reduced 
dimension  CSOM  matrix 


R  *  DHRD  (38) 

Is  defined.  For  both  the  beam  space  and  SA  space  processes,  the  secondary 
beam  output  variance 


2  -H  =  - 

ob  *  w  R  w  (39) 

Is  minimized  with  respect  to  either  the  (L+l)  or  (N/P)  dimensional  vector 
w.  The  constraints  that  the  element  in  location  ( L/2) -el  of  w  be  unity  for 
the  beamspace  and 


N  =  wV  d  (40) 

for  the  SA  space  procedures  are  required  to  satisfy  the  distortionless 
signal  constraint. 

It  Is  a  direct  procedure  to  obtain  the  following  AGI  expressions 


mi.  - 1  ♦  r-  l“M  , 

1  *  r[Ltl]l 


and 


r  1(1  -t) 
AGIs  - 1  *  —rFi; 


(41a) 


(41b) 
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for  the  beam  space  and  SA  counterparts,  respectively,  to  Eq.  (32)  which 
corresponds  to  the  fully  optimum  element  space  configuration.  In  Eq.  (Ala), 


j 


s 


(42) 


where  *k  *  Id^d^/N  2(0  <  lk  <  1)  is  the  relative  response  level  of  the 

Interference  In  the  k-th  beam  output  of  the  beam  space  patch.  The  quantity  I 
can  be  thought  of  as  the  interference  response  level  averaged  over  all  (Ltl) 
beams  in  the  beam  space  patch.  In  Eq.  (41b),  |$  (0  <  ls  <  1)  is  the 
relative  response  level  of  the  Interference  in  the  SA  beam  output.  Note 
that  in  the  limiting  case  for  SA  MVDR  processing,  P=1  corresponds  to  only 
one  sensor  per  SA.  Here- the  SA  has  an  omnidirectional  response  so  that 
=  1  and  the  result,  is  the  same  as  Eq.  (32)  as  expected.  It  is  observed  that 
for  the  two  suboptimum  MVDR  processes  to  perform  equally,  then 


J 


> 


=  [L+l )£ 


(43) 


and  optimality  is  approached  only  to  the  extent  that  e  approaches  unity. 
Figure  7  gives  the  AGI  metric 


AGI(e) 


rfiTe  -  t) 
’  1  +  re 


for  the  suboptimum,  reduced  dimension  beam  and  subarray  space  MVDR  processes 
for  several  values  of  e.  The  same  values  of  the  interference  response 
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AGI 


CBF  MRA  INTERFERENCE-TO-NOISE 
RATIO  (dB) 


Figure  7.  Comparison  of  Suboptimum  (Subarray  and  Beam  Space) 
and  Optimum  (Element  Space)  MVDR  Process  AGI 
Versus  Interf erence-to-Noise  Ratio. 
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level,  1,  for  the  full  aperture  CBF  are  used  as  In  Figure  5.  It  is 
significant  that  at  high  Interference-to-nolse  levels  (r)  the  suboptimum 
procedures  are  nearly  equivalent  to  the  optimum  element  based  process  except 
for  I  »  1/2.  Furthermore,  It  Is  primarily  only  for  large  r  that  substantial 
sldelobe  Interference  AGI  Is  obtained.  For  malnlobe  PIN,  when  1  *  1/2  the 
beamspace  MVOR  would  have  a  substantial  performance  loss.  This  Is  because  e 
only  differs  from  1/2  by  the  average  sldelobe  level  of  the  Interference  over 
the  remaining  beams  In  the  patch  and  this  would  be  a  small  number.  Thus, 
for  Interference  within  the  malnlobe  subarray  MVDR  would  be  superior. 


CONCLUSIONS 


The  fundamentals  of  adaptive  beamformer  (ABF)  Implementation  using 
systolic  computing  arrays  has  been  presented.  It  has  been  shown  that  for  a 
continuously  updating  ABF,  a  direct  open  loop  realization  can  be  obtained 
with  a  linear  systolic  array  consisting  of  just  two  types  of  functional 
computing  cells.  The  performance  of  a  generic  ABF  system  has  been 
reviewed.  Finally,  the  problem  of  computational  burden,  memory 
requirements,  and  extreme  convergence  time  associated  with  arrays  having 
large  numbers  of  elements  has  been  addressed  by  showing  that  systolic  array 
techniques  need  be  applied  only  at  the  second  stage  of  a  cascaded 
beamformer.  The  tremendous  saving  In  MVDR  Implementation  hardware  with 
application  of  the  suboptimum  processes  could  offset  the  loss  of 
performance.  The  preferred  suboptimum  MVDR  processes  use  subarray 
prebeamforming  because  it  Is  extremely  regular  in  its  architecture;  it  is 
not  frequency-dependent;  and  it  yields  better  A6I  performance  for  the  same 
MVDR  complexity. 
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OTIC  12 

APL/Johns  Hopkins  1 

APL/U.  Washington  1 

ARL/Penn  State  1 

ARL/U.  Texas  1 

MPL  Scrlpps  1 

Woods  Hole  Oceanographic  Institution  1 

ESL/TRW,  Inc.  (B.  Chaterjee)  4 

Contract  No.  N00001 40-83 -C -0110 
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